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Abstract — We consider the age-of-information in a multi-class 
MjG/l queueing system, where each class generates packets 
containing statns information. Age of information is a relatively 
new metric that measures the amount of time that elapsed 
between statns updates, thus accounting for both the queueing 
delay and the delay between packet generation. This gives rise 
to a tradeoff between frequency of status updates, and queueing 
delay. In this paper, we stndy this tradeoff in a system with 
heterogenous users modeled as a multi-class M/G/l queue. To 
this end, we derive the exact peak age-of-Information (PAoI) 
profile of the system, which measures the “freshness” of the status 
information. We then seek to optimize the age of information, 
by formnlating the problem nsing quasiconvex optimization, and 
obtain structnral properties of the optimal solution. 

I. Introduction 

Realtime status information is critical for optimal control 
in many networked systems, such as sensor networks used 
to monitor temperature or other physical phenomenon IT]; 
autonomous vehicle systems, where accurate position informa¬ 
tion is required to avoid collisions El; or wireless networks 
where realtime channel state information is needed to make 
scheduling decisions 0. In all of these systems, what matters 
is not how fast the update information gets delivered, but 
rather, how accurately the received information describes the 
physical phenomenon being observed. 

Recently, proposes the notion of age-of-information 
(Aol), which measures the average time between the gen¬ 
eration of an update message until it is received by the 
control unit; thus measuring the “freshness” of the available 
information. The early works on age-of-information consider 
homogeneous systems, where all entities use the same amount 
of resources for status update, and the length of status mes¬ 
sages can be modeled using an i.i.d. exponential distribution. 

In this paper, we consider a heterogenous systems where 
entities generate status messages with different length (service 
time) distributions. In particular, we consider a multi-class 
M/G/l queueing systems; where each entity generates status 
update messages according to a given distribution, and derive 
the exact peak age-of-information (PAoI) value for each entity, 
which is the average maximum elapsed time since the latest 
received update packet is generated, and captures the extent 
to which update information is delayed. We then consider a 
system where, for packet management, at most one packet can 
be kept in the system, and compute the PAoI for the multi¬ 
class M/G/1/1 queue. 


Next, we turn our attention to the problem of optimizing the 
PAoI by controlling the arrival rate of update messages (i.e., 
the sampling rate of the physical process being observed). We 
formulate the optimization problem as a minimization of a 
quasiconvex cost function of the system age-of-information 
prohle. We show that in the M/G/1/1 case, the optimization 
problem is a quasi-convex program and derive properties of 
the optimal solution. In the general M/G/l case, however, 
the problem is a general non-convex program for which we 
derive an approximate solution. 

The notion of age-of-information was first introduced in |4| 
in the context of a single source modeled as an M/M/1 queue, 
and extended to multiple sources in E). In 0, Q the authors 
consider the problem of minimizing the age of system state 
in vehicular networks. In fH the Aol is analyzed for a system 
with random delays, and in l9] PAoI is derived for a single¬ 
class M/M/1 queueing system. 

Our work differs from these prior works in a number of 
ways. First, we focus on the PAoI metric, which is closely 
related to the Aol, but is much more tractable, thus facilitating 
its optimization. Second, we consider general service time 
distributions, whereas previous works focus mainly on expo¬ 
nential service time. Finally, we minimize the cost of the PAoI 
in a system with heterogenous service requirements; where the 
service requirements of different entities are modeled using 
quasi-convex cost functions of the PAoI. 

It is important to note here that Aol is different from the 
traditional delay metric considered in communication systems. 
Indeed, our results show that PAoI minimization is equivalent 
to minimizing the sum of update interval and update packet 
delay. Due to this difference, the ultimate optimization prob¬ 
lem turns out to be non-convex. 

This paper is organized as follows. In Section]^ we present 
the system model. We derive the PAoI values for M/G/l 
and M/G/1/1 in Section III We consider the system cost 
optimization problem in Section IV and present numerical 


results in Section |V] We conclude the paper in Section VI 


II. System Model 

We consider a system that consists of a set of N entities, 
denoted hy Af = {1,...,A}. To disseminate entity status 
information, the system regulates how frequently each entity 
updates its status. We denote this decision by an update rate 
vector A = (Ai,..., \ n ), where A„ is the update rate of entity 





n. We assume that the update process is Poisson with rate A„. 
After generation, update packets are relayed to a central unit 
for processing. To ensure that each entity eventually updates 
its status, we require that 0 < Amin < A„ < Amax for all n. 
We denote by A = {A : Amin < An < AmaxjVn} the set of 
feasible rate vectors. 

A. Single queue model 

In a practical system, different update information streams 
will share the limited system resources for delivery. We model 
this by the system shown in Fig. where all update messages 
go through a single server queue. This single server queueing 
model is the same as that adopted by prior work on Aol H, 
0 , 0 . 



Fig. 1. The update packet delivery process in a 2-entity system. 

In this queueing system, each new class-n packet arrival 
to the queue represents the generation of a new entity n 
update packet. Hence, class-n packets arrive according to a 
Poisson process with rate An. Departures from the queue, 
on the other hand, represent update packet reception events 
at their destinations. To model the heterogeneous resource 
requirements, e.g., entities may have update messages with 
different length distributions, we assume that each update 
packet from entity n requires a random service time which 
is i.i.d. with mean Xn and second moment ?/„. For convenience, 
we denote Xmax = max„ Xn and j/max = max„ yn- 

This system is indeed equivalent to a multi-class MjGjX 
queue. This model captures two key features of communi¬ 
cation in networked systems: (i) resources are shared among 
different update streams, and (ii) queueing can occur during 
traffic delivery. Adopting this simple model allows us to focus 
on the age aspect of the update information. 

B. Age-of-Information 

The status age of an entity at a particular time instance t is 
dehned to be the time elapsed since its latest received update 
packet was generated. Fig. shows the status age (denoted by 
An{t)) of entity n. The dropping points of A„(f) are the time 
instances when an update packet is received, which resets the 
age value to a lower level (i.e., the current time minus the 
generation time of the new update packet). 



'n1 '^n2 


Fig. 2. The evolution of the status age of entity n in the system. Here 
denotes the k-th peak of age. 


Given A„(f), the average status age of entity n is defined 
as: 

1 A 

Aav,n = lim — / An{t)dt. (1) 

T^oo T 

This metric is called the age-of-information (Aol) and was 
hrst considered in HI. However, the Aol metric is hard to 
analyze. Moreover, in many systems, it is often the maximum 
status information delay that determines the performance loss 
0. Thus, we instead focus on the average peak status age. 
Specifically, let Ank denote the fc-th peak value of A„ (t) (See 
Fig. I). The peak age-of-information (PAoI) metric A„(A) is 
defined as: 

1 ^ 

^p,n(A) A lim yp'^Ank- (2) 

Here we explicitly express PAoI as a function of the update 
rate vector A. The PAoI metric was hrst considered in jJl 
for the M/M/1/1 queue. It represents the maximum age of 
information before a new update is received. PAoI is closely 
related to the previously considered Aol metric Aav,n, e.g., 
in a and 0, but is much more tractable, thus facilitating its 
optimization!^ 

C. Optimizing update rates 

We model the system performance using a function of 
PAoI, to capture the fact that delay in status update often 
causes a proportional performance degradation. Specihcally, 
we consider the following system cost function, i.e., 

C'sys(A(A)) A maxC'„(Ap„(A)). (3) 

n 

Here C'„(A„) is the cost of entity n for having a PAoI of A„, 
and it is assumed to be a quasiconvex and non-decreasing 
function in with C„(0) = 0. Our objective is to hnd 
an update rate vector A G A to minimize C'sys(A.(A)), i.e., 
minimize the maximum cost over all entities. 

III. Computing PAoI 

In this section, we compute the PAoI for two different 
queueing models, i.e., the M/G/1 model and the M/G/1/1 
model. The two models differ in the way update packets are 
managed. In the M/G/1 model, new packets are queued if 
the server is busy, while they are discarded in the M/G/1/1 
model. This packet management scheme was also considered 
in 0 for the M/M/1 model. 

A. A general result for GjGjl queues 

We first derive a useful result for general GjGjl queues. 
Denote by the inter-arrival time of entity n packets, and 
let Wn be the waiting time of entity n packets and let T„ be 
the total sojourn time in the queue. We have the following 
proposition, in which the superscript “ggl" is used to indicate 
the relationship to GjGjl queues. 

Proposition 1: In a GjGjl queue, the PAoI is given by: 

ASS/=E{J„ + T„} =E{/„ + A„ + IV„}, Vn. o (4) 

'intuitively, PAoI provides an approximate upper bound for Aol, as it only 
considers the average sampled at the peak moments, whereas Aol in ^ 
computes the time average value of instantaneous age. 













Proof: This relation follows from Fig. where we see 
that the PAoI is equal to the time from the generation of an 
update packet until the completion of the next update packet, 
plus their inter-arrival time. ■ 

Equation 0 shows that PAol is indeed the sum of the 
update interval and update packet delay. For a multi-class 
G/G/l queue, we also know that the Aol is given by ||4ll : 

= + (5) 

Comparing 0 and Q, we have: 


/isffi _ 499 ! 

-^av,n -^p,n 


( 6 ) 


= A„ + y)} - e{/„}e{j„ + r„} 


— A„ f E|/„r„} — Ej/njElTn} + 






Equation (j^ provides a way for checking how close the two 
metrics are to each other. It can be seen that when is a 


constant, i.e., periodic arrival, We will 

see in the next subsection that jn the single 

class M/M/1 case. Thus, PAoI serves as an upper bound for 
Aol. More generally, the following lemma shows that PAoI 
approximates Aol for general single-class G/G/l queues. 
Lemma 1: For a general single-class GjGjX queue, 

} _ ae{/}^ < Alf < Af 1 -f AE{J2}/2.0 
Proof: See Appendix A. ■ 

B. PAoI for multi-class MjG jX queue 

Using Q, we can compute the PAoI for each entity in the 
MjGjX queue. 

Proposition 2: The PAoI for a multi-class MjGjX system 
is given by: 


Amgl _ 

■^p,n \ 

'^n 




o 


(7) 


2(1 — ^jXj) 

In (0 we have used the P-K formula to compute the waiting 

time in the MjGjX queue ifTOl . It is necessary to ensure p = 

< 1 , so that the queue is stable and the PAoI is finite. 

We can use 0 to compute the PAoI for a single class 

M/M/1 system. In particular, with N = 1 and exponential 

service time of rate p, 0 becomes: 

1 / 1 P , 

^ ' ( 8 ) 


mml 


=-!+-+ 


A 

p py ' p ' I- 
In contrast, the Aol for M/M/1 derived in ID, is given by, 

A— 1 = + (9) 

P\ P ^-Pj 
Comparing 0 to Q, we have, 

A™-1-A-™i = -p=4- (10) 

P P 

Hence, PAoI is a close upper bound of Aol for M/M/1 
queues, yet is much more tractable. 

Conservation Laws for PAoI: From Q, we also obtain 
the following conservation formula for PAoI: 


^A„A™f = N + p + NW. 


( 11 ) 


We also see from Q that: 

'^n '^m 

which implies, 

4 - ^ - Xn) - {A^A _ ( 12 ) 

Hence, the relationship between A™A ^nd A^A completely 
determined by A„ and Am- For example. Am > A„ implies 

Amgl _ ^ Amgl _ 

p,n •^n p,m -^rn- 

C. PAoI for M/G/l/l queue 

Let us now consider the case when the server does not queue 
incoming update packets. This can be viewed as the server is 
performing packet management i). 

Proposition 3: The PAoI for a multi-class M/G/1/1 is 
given by: 

+ 4 + O (13) 

Proof: Let Zjn be the expected time to complete service 
for a class n packet starting from the moment a class j packet 
begins receiving service at the server. We have: 

+ (14) 

j 

Equation ( [T4| can be understood as follows. Since the status 
age of an entity decreases only when its next update packet 
is served, the (expected) peak age can be broken down into 
three components (see Fig. 0. The first component Xn in 
( pAl is the processing time of the current update packet. The 
second component is the expected time needed to get the 
next arrival (because packets arriving during busy periods are 
dropped). Then, the third component is the expected time 
needed until the completion of the next class n update packet 
(the third term), which is Zjn if the next arrival turns out to 
be a class j packet, resulting in an average time of ^Zjn- 



Time to serve 
the next packet 


Fig. 3. Evolution of entity n’s status age in the M/G/1/1 system. 

We now solve for A™A^. Note that for any i 7 ^ n, we have: 
1 ^ A,. 

(15) 


Zji -^j+ X 




A 


Thus, Zji — Zki = Xj — Xk for any j, k i. Plugging this into 
( [T5| ) and using the fact that Za = Xi, we obtain: 

+ - -f yXi + ^ ~ 


^ji — 


Therefore, 


1 


k^i 


Zji — Xj^ — P Xj ^ / — 

1 _L- ^ 

k^i 
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\ ^k^k 

r^' 


















Using this in (14 1 , we get: 

Amgll _ i 

■^p,n — ~ A ~ A 


j¥^n 


1 AfcXfc 

XT + 


A. 


— 2x7t, 4“ 


1 , Thk^i ^kXk 


\ Amgll _ \ 

^n-^p^n 


Xrt — 
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H” 1 ■ 


Similar to (121, the relationship of and A, 


mgll 


uniquely determined by A„ and A„ 
lowing conservation formula: 


p,n 

. We also have the fol 


E \ A^mgll _ 


N + {N+l)p. 


Comparing Q and ( [T3| ), we have: 

z- ^k^k 


Amgll _ Amgl _ 
■^p,n -^p,n 


E, Aj% 


Xr 


2(1 El XjXj) 


mm : 
A 

S.t. 


C^sys(A) — mux Cti (xyi A 

n 

A e A. 


A. 


The following lemma shows that although (20 1 is not convex, 
it can still be efficiently solved. 


Lemma 2: Problem ([20|) is a quasiconvex program. O 


Proof: First, we see that = Xr, 


1 + 


An An _ 

Rearranging the terms in the above gives ( [13| ). ■ 

It is interesting to note that © does not require the second 
moment of service time, which is generally required in the 
analysis of MjGjX queues. This can be attributed to the fact 
that in the MjGlXjX system, packets are never held in the 
buffer, thus the residual service time does not play a role in 
the computation of 0- Also, when = 1, 0 recovers the 
result from m for the MlAIjljl queue. It is also interesting 
to see in ( [T3| ) that due to packet discard, the constraint p < 1 
can actually be violated. 

We note from 0 that, for any achievable PAoI vector 

Amgll _/ A^ngll 


conclude that C'sys(A) is also quasiconvex in A and (20i is 
a quasiconvex program over the convex set A ifTTl . ■ 

Therefore, the optimization problem ( |20)i can be solved by 
the bisection procedure described below lUTl . Define 

0 C'sys(A) <t 

oo else. 

We see that C'sys(A) < f is equivalent to ft < 0, i.e., if A 
ensures ft < 0, it also ensures Csys(A) < t. Hence, we then 




use the following bi-section algorithm to solve (20i. 
Bisect ion: Set I = 0 and u = max„C„(a;i] 


Amin 

u — I < e: 


. Fix e > 0. Then, repeat the following until 


1) Set t = {I + u)/2 

2) Solve the following problem: 

min : 1, s.t. ft < 0, A S A. (22) 

3) If (|2^ is feasible, set u = P, otherwise set I — t. O 


(16) 


Using 0. we also obtain the following properties of the 
optimal solution A* to (20 1 , where E denotes entrywise larger. 


Since the right-hand-side (RHS) does not depend on n, this 
implies that: 

\ Amgll _ \ m _ \ Amgll _ \ „ /17', 


IS 


Lemma 3: Let A be an optimal solution of (20i. Then, 
(i) 3 A E A* such that max„ A„ = Amax and Csy^A) = 
C'sys(A*), and (ii) if all entities are identical, i.e., xi = X 2 
and Gn{A) = GmiA), then A„ = Amax, Vn is an optimal 
solution. O 


(18) 


.(19) 


Proof: First, note from (17 1 that for a given A, if we let 
no = argmin„(Aji — x„), then we can express each Am as: 

, _ , — Xuq 


v? 


(23) 


This shows that can be much smaller than when 

the update rates are large, i.e., when p is close to 1. Thus, even 
though packets can be dropped in the MjGjXjX system, such 
dropping may actually result in PAoI reduction as queueing 
delay is reduced. 

IV. PAoI OPTIMIZATION 

Having computed the PAoI for the two cases, we now con¬ 
sider the problem of optimizing the update rates, i.e., minimize 
C'sys(A). This formulation enables us to provide differentiated 
service to different applications. In the following, we start with 
the M/G/l/l queue and then consider the M/G/1 queue. 

A. MIGlljl optimization 

In this case, the utility optimization problem takes the 
following form: 

1 + Efc XkXk , 


° A -T 

Consider an optimal solution A* and let A* be the corre¬ 
sponding PAoI vector. We construct a A as follows. Denote 
tiq = argmin„(A* — Xn). Then, we keep the ratio between 
any pair of rates fixed and proportionally increase all A„ until 
An* = Amax- From (17 1 and (13 1 , we see that A < A*. Since 


each G„(A) is nondecreasing, we have Gn{An) < Gn(A*), 
which implies Gsys(A) = Gsys(A*). This proves (i). 

When all entities are identical, the rates must be the same 
for all entities. Hence, An = Amax for all n. ■ 


B. MIGjl optimization 

In AljGjX, the optimization problem becomes: 
min : G^y^ (A) = max Gn() 

1 , , Ej XjPj 

— -I- X.. -I- -- 


(24) 


s.t. 


2(1-E. V,™ 


-) ( 20 ) 




^3 -'3-^3 

< 1, An > 0. 


IS a 


Different from problem ( |20) i, here the LHS in constraint ( |25| ) 
is a sum of two linear-fractional functions, which may not 
be quasiconvex any more. Thus, to proceed, we approximate 
An (A) with another function Bn{\) defined as: 

1 Ej XjVj 


linear-fractional function in A. Since each Gn(An) function is 
quasiconvex and nondecreasing in An, Gn(An) is quasiconvex 
in A. As the max operator preserves quasiconvexity, we 


Bn{X) = 2 max ( --h x. 

An 


2(1 El AjXj) 


(26) 


That is, we solve problem (24 1 with A™A replaced by Bn- The 
main advantage of introducing Bn{X) is that it is quasiconvex 


























in A. Hence, the function C'sys(-B(A)) = max„C'„(B„) can 
be efficiently minimized by the Bisection algorithm. 

We now look at the performance of the approximation. 
Define /3„ the maximum increasing slope of Cn{A), i.e., 

/3„^inf{/?: |C„(Ai)-C„(A 2 )| <^|Ai-A 2 |, 

We then have the following lemma, where A^ is the optimal 
solution of the approximation program. 

Lemma 4: Let A* be an optimal solution of the original 
problem ( |2^ and denote A* the resulting PAoI. Then, 

Csys(A*) < CsysCA^) < + max/3„A^ O (27) 

n 

Proof: See Appendix B. ■ 

When each C'„(A) is linear in A, i.e., C'„(A) = WnA, we 
have Pn = Wn- In this case, we can conclude from ( |3^ that 

Csys(A*) < C,y.{X%) < 2C,y,{X*). 

V. Numerical Results 

We present a simple numerical example with N = 2 entities, 
and constant service times Xy = 1 and X 2 = 3. Amin = 0.01 
and Aniax = 10. The cost functions are given by C'i(Ai) = 
4Ay and < 72 (^ 2 ) = A 2 , and C'sys(A) = max(Ci,C 2 ). 

In Fig. 1^ we plot the cost values for the MjGjlll queue. 
The minimum value of Csys(A) is achieved at Ai = 10 and 
A 2 = 6, with a resulting Cy = 60.84 and C'sys(A) = C 2 = 
61.36. In this case, the PAoI vector is A = (3.9, 7.83) and the 
results match Lemma 0 



Eig. 4. Cost value for the M/G/1/1 queue. 


We next look at the M/G/l case in Fig. Since p < 1 
must be satisfied to ensure a finite PAoI, if a A violates p < 
1, we set its PAoI value to a constant (the flat region). In 
this case, the minimum is achieved at A = (0.29,0.125) with 
A = (6.56,13.11). Thus, G 2 = 171.92 and Cy = C,ys(X) = 
172.15. It can be verified that ( [T^ holds. 



Eig. 5. Cost value for the MjGIl queue. 


We also compute the optimal solution of the approximation 
approach to be Ag = (0.285,0.17). The resulting PAoI vector 
is A = (8.94,13.31) and the cost vector is (^ 1 ,^ 2 ) = 
(319.69,177.16), implying Gsys(A^) = 319.69. It is not hard 
to verify that Gsys(A^) < 2Gsys(A*). 


VI. Conclusion 

We study the age-of-information in a general multi-class 
M/G/l queueing system. The age-of-information is a new 
metric for system performance that represents not just the 
queueing delay, but also the delay in generating new infor¬ 
mation updates. Our main contribution is to generalize the 
available results to systems with heterogeneous service time 
distributions, accounting for the fact that different entities may 
have different service requirements for their status updates. 
We derive exact peak-age-of-information expressions for both 
a M/G/l system and M/G/1/1 system. Using the PAoI 
measure allows us to optimize system cost, as a function of 
PAoI, by choice of the update interval. 
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Appendix A - Prooe oe Lemma[T] 

Here we prove Lemma [T] We drop all subscripts as W = 1. 
Proof: Since T = W + X and X is independent of I, we 
have E{/„r„} - E{/„}E{r„} = E{lW} - E{/}E{1U}. 
Denote Tq the time it takes to clear the packets in the queue 
when a new packet arrives (not including the new arrival). 
Then, W = {Tq — /)'*'. Here I is the inter-arrival time until 
the next packet arrives. Since I is independent of Tq, we have; 

E{IW} 

= [ r i{t-1) f{t)dtf{i)di 

JI Jt^i 







= 1(^1 ^ > 1}^ f{l)dl (28) 

< [if tf{t)dtf{l)dl. (29) 

Jl Jt^I 

Using (Tg — J)+ 1 >Tq and / > 0, we have: 

poo 

/ tf{t)dt<E{TQ}<E{{TQ-I) + }+E{l}. (30) 

Jt=i 

Plugging this into ( |29] l, we get; 

E{/M^} < j^l{E{iTQ-!)+}+E{l})f{I)dI 

= E{/}E{(rQ-/)+}+E{/}^ (31) 

Using this in (j^, we obtain: 

To derive the lower bound, we have from ( |28| l that: 
E{/1U} = ^ (^I J^°[^tf{t)dt-I^PT{TQ >I}^f{I)dI 

> I (^E{Tq} - tfit)dt^ f{l)dl - E{/2} 

> /(^E{Tq} - f{l)dl - E{/2} 

> E{J}E{VU}-2E{J2}. (32) 

In the last inequality, we have used Tq > W. Plugging (32i 
into '0 proves the lower bound and completes the proof of 
the lemma. ■ 

Appendix B - Proof of Lemma|4] 

We prove Lemma Inhere. 

Proof: From the dehnition of Bn, we see that given any 
A, < BniX) < 2A™si(A). Thus, we have for 

that: 

ay,{Ai\*B)) < ay^iBiXD) (33) 

< ay.(B(A*)) < a,,(2A(A*)). 

Using the dehnition of we have for each n that: 

Cni2Ani\*)) < Cn{An{X*) + PnA^. (34) 

Taking the max over the above and combining it with (33i, 
we see that (|27ll follows. ■ 





